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Quark and gluon form factors to three loops 
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We compute the form factors of the photon-quark-anti-quark vertex and the effective vertex of 
a Ifiggs boson and two gluons to three-loop order within massless perturbative Quantum Chromo- 
dynamics. These resuhs provide building blocks for many third-order cross sections. Furthermore, 
this is the first calculation of complete three-loop vertex corrections. 
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In the recent years various next-to-next-to-leading or- 
der (NNLO) calculations to physical observables have 
been completed. Among them are the total threshold 
cross section for top quark pair production in electron 
positron annihilation the Higgs boson production 
in gluon fusion 0, H, the rare decay rate of the B 
meson into a meson containing a strange quark and a 
photon (sL [yl and the three-jet cross section at lepton 
colliders [3, @]- There exist also a few results at next- 
to-next-to-next-to- leading order (NNNLO), like the to- 
tal hadronic cross section in electron positron annihila- 
tion J9|, the hadronic t lepton 9] and Higgs boson de- 
cay It is common to all NNNLO results that the 
calculation can be reduced to two-point functions and 
that only one mass scale is involved in the computation. 

In this Letter we provide the first NNNLO calculation 
of a three-point function within Quantum Chromody- 
namics (QCD). To be precise, we consider gauge invari- 
ant building blocks for NNNLO cross sections, namely 
the virtual third-order corrections for the hadronic Higgs 
boson production and the process e^e^ — > 2 jets. The re- 
sults are conveniently expressed in terms of form factors 
of the photon-quark and the effective gluon-Higgs boson 
vertex originating from integrating out the heavy top- 
quark loops. Denoting the corresponding vertex func- 
tions by and F^" , respectively, the scalar form factors 
are obtained via 



{qi ■ q2 9^iu - ql,^i q2., 



2(1-6) 



g ■ 



(1) 



where d = 4 — 2e is the space-time dimension, q = qi + q2 
and qi (52) is the incoming (anti-)quark momentum in 
the case of Fq, and Fg depends on the gluon momenta 
qi and 52 with polarization vectors s^{qi) and £'^(52)- 
Some sample Feynman diagrams contributing to Fq and 
Fg are shown in Fig. [TJ Starting from three-loop order 
a new class of diagrams occurs, the so-called singlet dia- 
grams, where the external photon is not connected to the 
fermion line involving the final-state quarks (see Fig. [T] 
(b)). Since at three- loop order there are no counterterm 
contributions to the singlet diagrams and furthermore 
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FIG. 1: Sample Feynman diagrams contributing to the Fq 
((a) and (b)) and Fg (c) at three-loop order. Straight and 
curly lines denote quarks and gluons, respectively. 



there is no corresponding real emission contribution the 
sum of all diagrams has to be finite. This constitutes an 
important check on the correctness of our result. 

In the recent years the evaluation of the three-loop 
form factor has attracted much attention. After the pio- 
neering work more than 20 years ago [111, [13, [3 where 
the quark form factor has been computed to two-loop or- 
der the corresponding quantity for the Higgs-gluon cou- 
pling has been evaluated by Harlander in Ref. (see 
also [15]). The latter constitutes a building block for 
the NNLO predictions of the Higgs boson production in 
gluon fusion at the Fermilab Tevatron and CERN Large 
Hadron CoUider [1,[1,I3. More recently, in Ref. ^ the 
two-loop results have been reconsidered and more terms 
in the e-expansion have been added in order to match 
the three-loop accuracy. Furthermore, in Refs. [l3, EB] 
almost all master integrals necessary for the three-loop 
calculation have been evaluated. However, the most com- 
plicated master integrals are still unknown. 

First steps towards three-loop results for the form fac- 
tors have been undertaken in the Refs. [l^ where 
the pole parts of Fq [l^ and Fg have been extracted 
from the behaviour of the three-loop coefficient function 
for inclusive deep-inelastic scattering [21.] . Furthermore, 
in Ref. 20] also the finite part of the fermionic contri- 
bution to Fq could be evaluated. With our calculation 
we were able to confirm these results but also add the 
finite contributions which are necessary for the physical 
observables. 
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FIG. 2: Three most complicated master integrals entering the 
result for the three-loop form factor. The notation is adopted 
from Ref. [13, El. 
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For the evaluation of the Feynman integrals we de- 
veloped two independent set-ups which have in common 
that a reduction of all occurring integrals to so-called 
master integrals is performed in d space-time dimen- 
sions. Afterwards the (e-expanded) master integrals are 
inserted. 



Following Refs. [22, l23|, |2j| one considers integral rep- 
resentations of the coefficient functions of the individ- 
ual master integrals in the limit of large space-time di- 
mension d, evaluates several expansion terms and re- 
constructs in this way the complete rational dependence 
on d. The most CPU-consuming step, the large d ex- 
pansion, has been performed by a program written in 
ParFORM [HI, [2^, the parallel version of the computer 
algebra program FORM [27'|. For the singlet contribution, 
which involves the most complicated integrals, also a sec- 
ond approach has been employed. After generating the 
Feynman diagrams with the help of QGRAF [l^ they are 
further processed with q2e and exp 29, 3^ where a map- 
ping to the underlying family of the diagrams is achieved. 
In a next step the reduction of the integrals is performed 
with the program package FIRE [Slj which implements 
a combination of the Laporta algorithm and a gen- 
eralization [3^ of the Buchberger algorithm (see, e.g., 
Ref. [13]) to construct Grobner bases. 



Our results are expressed in terms of 22 master in- 
tegrals. Eight master integrals are either two-point 
functions or products of one- and two-loop integrals 
and are thus well-known since many years (see, e.g., 
Ref. [il,!!!!!!,!!^). The results for eleven three-point 
master integrals can be found in Ref. [13, EBl, however, 
the three most complicated integrals, which are shown in 
Fig. [21 are not yet known in the literature. Our calcu- 
lation fixes, by comparing with Ref. [l9j . the divergent 
parts of Aq^2 and Ag^^ and the finite part of Ag i and 
leaves only three coefficients of the e expansion unde- 
termined. The results read (assuming massless propa- 
gators of the form 1 / {k^ + iO) and pulling out a factor 
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We obtained a numerical result for the coefficient Xg^i 
using the MeUin-Barnes (MB) method [H, [H, HO] , start- 
ing from the general MB representation for the tennis 
court diagram of Ref. [4l[ , and applying the correspond- 
ing packages [H, H^. To evaluate numerically Xg^2 and 
X<)^4 we used the program FIESTA [H] which is a conve- 
nient and efficient implementation of the sector decom- 
position algorithm. Our results read 

Xg,! « 1429(1) ,Xg,2 « 528.0(4) , Xg,4 « -2085(5) , (5) 

where the accuracy is sufficient for all foreseeable physical 
applications. Finally, let us mention that we evaluate the 
colour factors with the help of the program color [45j . 

In the following we want to present explicit results for 
Fq and Fg. We parameterize the results in terms of the 
bare coupling which allows us to factorize all occurring 
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logarithms of the form ln((5^//i^) where = —(f > 0. 
Furthermore, we cast the results in the form (a; = q, g) 

n \ ^ / 

(3) 

and split Fq into the singlet, fermionic and remaining 
gluonic part 

i^f ) = i^f + Ff + ^ g,,i^f , (7) 

where n/ stands for the number of active quarks. The 
results for Fq^'' and Fg^-* (expanded in e sufficient for 
the three-loop calculation) can be found in Eqs. (3.5) 



and (3.6) of Ref. and F^^\ ^ and Fj^^^'"^ are given 
in Eqs. (7), (8) and (6) of Ref. [20|, respectively. The pole 
parts of F,^^^'^ and F^ ^ are listed in Eqs. (3.7) of Ref. [H 
and (9) of Ref. [l^l , respectively. Our expressions agree 
with all these results which constitutes a strong cross 
check since in Refs. [l^ [1^ a completely different ap- 
proach has been chosen to evaluate the Feynman inte- 
grals. In particular, no reduction to master integrals has 
been performed. In this Letter new results for F^f^'^, 

Ff and Ff ) are presented. Since the pole parts 
are already available in the literature we display only the 
corresponding finite parts which read in the case of a 
SU{Nc) colour group 
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where Cp = {N^ - l)/{2Nc), Ca ^ N^, T ^ 1/2 and 
^abc^abc. ^ (-^2 _ ;^)(^2 _ 4)/iVc. Inserting numerical 



values leads to Fq'^^'^^"^ Isr 



-13656. 
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and F, 



(3): 



26102.7 - 8298. 8n/ + 585. Sn^ ± 27Mg.i ± 

c one sigma 



21.6(59_2, where 5g,i — 1 corresponds to t^ii 
uncertainty given in Eq. ([5]). 

It is interesting to specify our result to a supersymmet- 
ric Yang-Mills theory containing a bosonic and fermionic 
degree of freedom in the same colour representation. This 
is achieved by setting Ca ^ Cf = 2T and n/ = 1 which 
leads to 
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Although we do not know three coefficients analytically, 
we believe that the growth of the transcendentality level 
continues when going to the next order in e so that all 
the results are at most of transcendentality six as was 
predicted in Refs. [U |4^ . It is interesting to note that 
these terms agree between the two form factors. 

To summarize, in this Letter we compute the form fac- 
tors of the photon-quark and effective Higgs boson-gluon 
vertex to three-loop order within massless QCD. Our re- 
sults constitute important building blocks for a number of 
physical applications. Among them are the two-jet cross 
section in e~^e~ collisions, the Higgs boson production 
in gluon fusion and the lepton pair production in pro- 
ton collisions via the Drell-Yan mechanism. Let us stress 
that our result represents the first complete evaluation 
of three-loop QCD corrections to a three-point function. 
Our results for the coefficients of the three master inte- 
grals 1 , ^9^2 and Aq^4 partially overlap with those of 
Ref. [47| where these integrals were evaluated in a direct 
way. Agreement has been found for all common coeffi- 
cients. 
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